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Using nonperturbative results obtained recently for an uniformly accelerated Unruh-DeWitt detec- 
tor, we discover new features in the dynamical evolution of the detector's internal degree of freedom, 
and identified the Unruh effect derived originally from time-dependent perturbation theory as oper- 
ative in the ultra-weak coupling and ultra-high acceleration limits. The mutual interaction between 
the detector and the field engenders entanglement between them, and tracing out the field leads to a 
mixed state of the detector even for a detector at rest in Minkowski vacuum. Our findings based on 
this exact solution shows clearly the differences from the ordinary result where the quantum field's 
backreaction is ignored in that the detector no longer behaves like a perfect thermometer. From a 
calculation of the evolution of the reduced density matrix of the detector, we find that the transition 
probability from the initial ground state over an infinitely long duration of interaction derived from 
time-dependent perturbation theory is existent in the exact solution only in transient under special 
limiting conditions corresponding to the Markovian regime. Furthermore, the detector at late times 
never sees an exact Boltzmann distribution over the energy eigenstates of the free detector, thus in 
the non-Markovian regime covering a wider range of parameters the Unruh temperature cannot be 
identified inside the detector. 

PACS numbers: 04.62,+v, 42.50.Lc 04.70.Dy 



I. INTRODUCTION 



The Unruh effect states that an observer while undergoing uniform acceleration in the Minkowski vacuum feels 
as if it lives in a thermal state at the Unruh temperature. It is usually demonstrated by way of time-dependent 
perturbation theory (TDPT) for a "particle detector" over an infinitely long duration of interaction [l], 0, H[ . Consider 
the detector as a quantum mechanical object with internal degree of freedom Q coupling to a quantum field $ by the 
interaction Hamiltonian 

Hj = AoQ(r)*(^(r)) (1) 

where Ao is the coupling constant, r is the proper time for the detector and z p (t) is the trajectory of the uniformly 
accelerated detector (UAD) with proper acceleration a. Suppose initially the detector-field system is in a product 
state, i.e., they are uncorrelated, 

| t — > -oo) = | E ) ® I M ) , (2) 

where | Eq) is the ground state of the free detector and | 0m) is the Minkowski vacuum of the free field. From TDPT, 
the transition probability Po-^ n from the ground state to the n-th excited state | E n ) of the detector is given to first 
order in 7 = Ao/8-7rmo, by [3| 
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which is nonzero. Here 

K(k) = J~ drdr'e-^-r ') ( o M |$( z (r))$(z(r'))|0 M ) = A_^_ (4 ) 

is the response function with r/ = J^^dr being the duration of interaction [20]. In particular, for a harmonic 
oscillator with renormalized natural frequency fl r , all the transition probabilities with n > 1 are 0{^ 2 ), so the only 
non-vanishing P of 0(7) is 

47rm e 2 ^-> - 1 1 ' 

From the Planck factor in (J5]) one can identify the Unruh temperature Tjj = ha/2irkB- When a — > 0, the transition 
probability per unit time P^i/rj vanishes, which implies that there is no excitation in an inertial detector initially 
prepared in its ground state Q. 

According to the above calculation of the transition probability one may expect that after the coupling is switched 
on, the reduced density matrix of the detector would evolve continuously from the ground state to higher excited 
states, with the detector finally ending up in equilibrium with a Boltzmann distribution ~ exp[— E n / (IcbTxj)] at the 
Unruh temperature Tjj. Then one can read off the Unruh temperature from the late-time distribution inside the 
detector. 

In this note we use the exact solutions to an Unruh-Dewitt detector 0, 0] in uniform acceleration interacting 
with a quantum field to show that the above described conventional scenario is good only in the Markovian regime 
corresponding to the limits of ultra-high acceleration or ultra-weak coupling. The transition probability calculated 
from the infinite-time TDPT is valid only in transient under restricted conditions. The evolution of the detector under 
general conditions is quite different from the above picture. 

The difference between this new understanding and the conventional picture originates from the fact that the full 
interplay between the detector and the field is included here whereas the conventional approach based on perturbation 
theory operative for infinite time ignores the backreactions - the use of TDPT over indefinite time amounts to invoking 
a Markovian approximation which imposes rather severe limitations. The conventional wisdom is built upon conditions 
which cannot reflect the most general features in the full dynamics of the detector-field system made possible here by 
the nonperturbative solutions. 

This paper is organized as follows: In Sec. II we consider a moving Unruh-DeWitt detector in 4-D Minkowski space 
coupled to a quantum scalar field $ and calculate the reduced density matrix after tracing out the field. In Sec. Ill 
we calculate the transition probability of the detector from the ground state to an excited state. We elaborate on the 
physical meaning of the two constants Ao, Ai in the theory. In Sec. IV we calculate the purity and the von Neumann 
entropy or the entropy of entanglement of the detector-field system and identify an effective temperature T e ff which 
is ostensibly different from the Unruh temperature. We identify the range of validity of the conventional results to 
that which corresponds to making a Markovian approximation with its strong limitations. In Sec.V we summarize 
our main results with remarks pertaining to the issues raised in this paper and conclude with a comment on how our 
nonperturbative results may bear on a new understanding of the Hawking effect of black hole radiance. 



II. REDUCED DENSITY MATRIX FOR DETECTOR 



Consider an Unruh-DeWitt (UD) detector moving in (3+1) dimensional Minkowski space. The total action is given 

by! 



(d T Q) 2 -n 2 Q 2 - I d i x^d^d ti ^> + A / dr / d 4 xQ(T)<P(x)6 4 (x» - z"(r)) , (6) 



where Q is the internal degree of freedom of the detector, assumed to be a harmonic oscillator with bare mass mo 
and bare natural frequency ^o. The scalar field $ is assumed to be massless, and Ao is the coupling constant. 

This UD detector theory behaves like the quantum Brownian motion (QBM) of a harmonic oscillator interacting 
with an Ohmic bath provided by the 4-D scalar quantum field ||. The QBM model is a useful comparison because it 
shows clearly the dissipative and stochastic behavior of the dynamics arising from the interplay between the system 
and the environment, and the influence functional treatment incorporates the backreaction of the environment on 
the system (which could be either the quantum field or the harmonic oscillator depending on what one is after) in a 
self-consistent way. Since the QBM model indicates that there are nontrivial activities at zero temperature [H, EI j 
we caution that even for the a = case the detector is not just laying idle but has interesting physical features due 
to its interaction with the vacuum fluctuations in the quantum field. 
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Unruh and Zurek Q have studied an exactly solvable QBM model where a harmonic oscillator interacts with a 
massless scalar field in 2-D. They derived the exact master equation for the reduced density matrix of the system 
(oscillator) at a temperature determined by the initial state of the field, and observed some general features different 
from the conventional Markovian results (valid for an ohmic bath at ultra- high temperature). One feature is the 
dependence of the UV cut-off in the master equation and the reduced density matrix, and thus also in the von 
Neumann entropy of the system. When the interaction is switched on, the factors in the master equation and 
the entropy have initial "jolts" over a very short time scale corresponding to the UV cut-off frequency, which cause 
significant change in the coherence of the quantum state of the system, while the coherence residing in each subsystem 
is rapidly transferred into the correlation between them. This effect is discussed in detail in [6J. Below we will see a 
similar behavior occurring in UD detector theory beyond the Markovian regime. 

Suppose the initial state of the system at r is given by a direct product of the ground state for Q and the 
Minkowski vacuum for $. In the Schrodinger representation, this initial state is a product of Gaussian functions, 



~AQ 2 - 



d 3 k d 3 k' 
(2tt) 3 (2tt) 3 



£ kk '$k$k' 



Owing to the linearity of the model, this state must evolve into a general Gaussian form 

Mr) = jV(r)exp \-A(t)Q 2 - B kk '(r)$ k $k' - 2G k (r)$ k Q 



(7) 



(8) 



after the coupling is switched on. Here we have used the DeWitt notation where each upper-lower indices pairing 
indicates an integration (or summation). So the density matrix for this pure state in the (Q, $k) representation can 
be written as p[Q, <&k, Q', ^ T \ = V'otQ, •I'k; t]i/jq [Q' , $' k ; r], and the reduced density matrix defined by tracing out 
the field $k reads 



J (Q,Q';r) = J V^ [Q,^;r]^[Q',^ k ;T] = exp [-G« (T)Q t Q 3 - F(t)] , 



(9) 



where i,j = 1,2, Qi = (Q,Q'), The factors G u and F could be obtained by solving the master equations @|. But 
here we do not need to solve them directly, because the two-point functions of the detector have been obtained in 
the Heisenberg picture in Ref.[3]. (These results are placed in Appendix A for convenience.) One can reconstruct 
the complete evolution of the reduced density matrix using those two-point functions of the detector by solving the 
simple algebraic relations 



G 11 
G 11 



G 22 + 2G 12 
G 22 - 2G 12 
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(P 2 )(Q 2 )-(P,QY 



(10) 
(11) 
(12) 



where (P,Q) = |( (PQ + QP) ) = (m /2)(d/dr) ( Q 2 ) (and ( [P,Q] ) = -ih). Substituting (IA1IIA4I) gives the 
values of G ZJ at every moment, and the factor F in ([9]) will be determined by a normalization condition. 

To compare with the transition probability the reduced density matrix ([9|) has to be further transformed to 
the representation in the basis of energy-eigenstate for the "free" harmonic oscillator Q: 



p R (Q,Q')= Yl Pm,n<t>m(Q)4>„(Q') 



(13) 



m,n>0 



where the wave function for the n-ih excited state is 

MQ) -- 



a 



2™n! A /7r 

with the Hermite polynomial H n (x) and the real constant 



H n (aQ)e 



-a 2 Q 2 /2 



(14) 
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Here f2 r is the renormalized natural frequency Q ■ Following the method given by Ruiz @, the matrix elements p„: „ 
could be extracted by comparing the coefficients of s™sV; on both sides of the following equation, 
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III. TRANSITION PROBABILITY 



The transition probability from the initial ground state to the first excited state is the m = n = 1 component of 
the reduced density matrix in energy eigenstate representation (|32l) : 



(P 2 )(Q 2 )-(P,Q) 2 ~(h 2 /A) 



3/2 ' 



{ [h-'a- 2 ( pa > + (ft/2)] [fto* ( Q3 ) + (k/2)} -(P,Q > 2 | 

Expanding the two-point functions of the detector (|AltiA4|) in terms of 7, ( Q 2 ) and ( P 2 ) look like 

ft 



Q 



pz 



2mof2 r 



-7<Q 2 ) (1) +0( 7 2 ), 



-m ft r + 7 ( P 2 



0(7 2 ), 



and ( P, Q ) - 0(7), which yield 
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(p 2 ) {1) + m 2 n 2 r (Q 2 ) (1) +o{Y). 



2hmoVl r 

When rj = t — To^>a _1 the approximate value up to the first order of 7 = Ao/87rmo becomes 



-1 \2 
R I V^a A Q 



'/ 



Ai + A - 21n(a/fi r ) 
2irQ r 



(20) 



(21) 



(22) 



(23) 



from (|AlllA4j) . Here Ao and Ai are large constants in (|A3j) and (|A4j) : Ai denotes the time resolution/frequency cut-off 
of this detector theory, while Aq denotes the time scale of switching on the interaction. 



A. Range of validity in conventional results 



We see that the first term of (|23p gives the conventional transition probability ([5]) from TDPT over infinite time. 
Only when Vt r r] ^> Ai, Ao, or a is extremely large, can the second term in ([23]) be neglected and ([5]) recovered. Hence 
the conventional transition probability §5§ is valid only in the limits of (a) ultra-high acceleration (a 3> f2 r and 
Ai <C arj <C a7 _1 ) or (b) ultra-weak coupling (a -1 , x Ai <C 77 <C 7 -1 )- 

Note that, in obtaining (|23p. we have assumed a -1 <C T) <C 7 _1 , when the system is still in transient. If a < 7, 
the conventional transition probability ([5]) has no chance to dominate at all. Mathematically, the first term in (|23|) is 
contributed by the poles at ±f2 + 17 or ±f2 — 17 in the K-integrations of two-point functions (see Eqs.(60) and (67) in 
[H). But when a < 7 the poles at ±ia would be closer to the real k axis than the poles at ±f2±ry, while the poles on 
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the imaginary axis k = ±ina, n € N become very dense for small a, so their contributions dominate the result and 
the first term of (|23[) becomes unimportant. 

In particular, the a = case is beyond the reach of TDPT over infinite time shown in Section I, and the conventional 
wisdom from perturbation theory that no transition occurs in an inertial detector is untenable. In contrast, our 
result indicates that the evolution of p^ 1 with a — behaves qualitatively similar to those cases with nonzero 
acceleration 0] . This agrees with the expectation from the observation that the UD detector theory is a special case 
of quantum Brownian motion [f| , where there is nontrivial interplay between the oscillator and the quantum field at 
zero temperature. 



B. Essential roles of Aq and Ai 



We see the presence of two constants Ao,Ai and may wonder whether they have some real physical meaning or 
are just part of a calculational tool or artifact. We will address these concerns here. To begin with, the presence 
of large constants corresponding to frequency cut-offs in the coincidence limit of two-point functions of the detector 
is a common feature for detector-field theories and quantum Brownian motion. For example, the Raine-Sciama- 
Grove(RSG) model in (1+1)D [l8|, in which the detector acts like a harmonic oscillator in a sub-Ohmic bath, also 
has two-point functions of the detector with dependence of large constants due to the infrared cut-off. 

Let us explore the physical meaning of Ao and Ai. Since Ao corresponds to the time scale of switching on the 
interaction, it could be finite in real processes, and for every finite value of Ao, the Ao terms in all two-point functions 
vanish at late times (see (|A3|) and (|A4[) 'l. Hence Ao will not be present in the late-time results. On the other hand, 
Ai is a constant of time, appearing from the very beginning and never decays. One way to see that it is a quantity 
of real physical meaning is that if Ai was subtracted naively, the uncertainty principle will be violated, namely, 
APAQ = -j\ P 2 ) ( Q 2 ) < Ti/2 at late times for a is small enough, as shown in FIG[T] 

Actually (f2"2"| is formally identical to the first-order transition probability from TDPT for a UAD with finite duration 
of interaction (to,t) [l(|, so is the result (f2"5)) . To verify this, first note that the 0(7) term in (|22[) is the 0(A§) term 
of 

jL-(E o ,0 M \Hq(v)\ Eo,0 M ). (24) 



Here Hq(j]) is defined by 



- , . 4- { T dsH(s) - -i f T dsH(s) , N 

H Q {r,) = e n ^o >H Q e (25) 



where H(t) is the total Hamiltonian for the combined system, Hq is the Hamiltonian for the "free" detector so that 
I Eq) is an eigenstate of Hq. Then in the interaction picture it is straightforward to show that 
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Pt 1 « / dn I dT 2 {E x - E ) (Ei |0(0)| Eo) e -*^-^Hn-«j ( Q M mz(n))<Z>(z(T 2 ))\0 M ) , (26) 

""r J to J to 

from which one recovers the finite-time transition probability from TDPT (cf. Eq.([3])) since E\ — Eg = h£l r - Math- 
ematically the large constants Ao and Ai are formally the same as the divergences found in Ref. [lOj. In Ref. 
it has been shown that these divergences can be tamed if one switches on and off the interaction smoothly, so can 
Ao. Nevertheless, the physical meaning of Ai here is totally different from those divergences of the finite-time UAD. 
Here we are looking at the real-time causal evolution problem ( "in-in" formulation) rather than a scattering transition 
amplitude ("in-out" formulation) problem. Also in our set-up we never turn off the coupling and Ai is present at 
every moment. 

As was already found in Ref. [7[ the constant Ai changes the scenario of the evolution of the system in an essential 
way. Ai is present whenever the coupling is on, so right after the initial moment the values of p^ n , m, n > jump 
from zero to large numbers depending on Ai, in a time scale indistinguishable from zero at the level of precision of 
this theory. This means that the initial distribution of pf n n peaked at the element Pq would upon the switch-on of 
the coupling collapse rapidly (rather than smoothly diffuse) into a distribution widely spread (if Ai is large) over the 
whole density matrix. After that the density matrix begins to redistribute itself and finally settles into a steady state. 
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IV. PURITY, ENTROPY AND EFFECTIVE TEMPERATURE IN DETECTOR 



After substituting the late-time two-point functions (|A7|1 and (|A8|) with Ai kept in ( P 2 ) into (fTB)) . we find that, in 
steady state, the diagonal terms of n in the free eigenstate representation do not assume the form of a Boltzmann 
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FIG. 1: Plot of (APAQ)/(h/2) at late times with Ai = 0. One can see that there exists a range of a (less than about 4) where 
APAQ = yj ( P 2 ) ( Q 2 ) is smaller than h/2, violating the uncertainty principle. So Ai cannot be naively subtracted. Other 
parameters in this plot are taken to be 7 = 0.1, fl = 2.3, mo = 1, and % = c = 1. 



distribution, while the off-diagonal terms of p^ n do not vanish at late times. Therefore the free eigenstates are no 
longer orthogonal when interactions are introduced. 

To define an effective temperature and the entropy of entanglement, one should choose another set of orthogonal 
states which have a similar form as the energy eigenstate of the free detector but with a in (TT5|) replaced by 

1/4 



P Z 



(27) 



Since our initial state is Gaussian, one can further apply a canonical transformation introduced by Unruh and Zurck 
to diagonize the reduced density matrix into 



~P R 



Hdr, 



{U - h/2) n 
' (W + h/2) 



n+1 ' 



where 



U 



yJ(P 2 )(Q 2 )-(P,Qf 



(28) 



(29) 



is the uncertainty function and U > h/2 is the Robertson-Schrodinger uncertainty relation. From (|28j) one immediately 
obtains the purity V for the reduced density matrix ([9]), 



V = Tr[p H ~p H ] 



h/2 
~TT' 



(30) 



In this case V 1 is the ratio of the uncertainty function U to its minimal value h/2 : hence one always has V € (0, 1) 
because of the uncertainty principle. Also the von Neumann entropy or the entropy of entanglement reads 



U 1 



U 1 



U 1 



U 1 



which is a measure of the entanglement between the detector and the quantum field. 
Further, in steady state, ( P, Q ) — * 0, so (|28l) becomes 



(31) 



Prn.n * """n 



y/(P*)(Q*)+ h/2 



n+1 



= x n R x p ~nhn r /k B T e tt 

— u mnPo A e 



(32) 



Now the diagonal terms assume a Boltzmann distribution with the effective temperature 



T, 



cif 



hQ r 



111 



y/(P 2 )(Q 2 )+ h/2 
V< P 2 )(Q 2 )-?i/2 / 



(33) 



This is not surprising: Even for such a simple system containing only two coupled harmonic oscillators, the ground 
state of the total system also looks like a thermal state with some effective temperature in view of the reduced density 
matrix for one of the oscillator 1121 . 
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FIG. 2: The evolution of the entropy of entanglement S(r/) in (|3ip for a detector initially in the ground state. At late times, 
the entropy converges to a large number if Ai is large. Here the parameters are taken to be a — 2, 7 = 0.1, £1 = 2.3, mo = 1, 
Ai = A = 10000, and h = c = 1. 



A. Realistic cases: fiAi ^> a, 7 



Since Ai corresponds to the cut-off frequency of the theory, in real processes it is most likely that ClA% 3> a, 7. 
The evolution of the entropy of entanglement in this case is shown in FIG|2] At the initial moment the entropy of 
entanglement has a sudden jump from zero to a large number ~ 0(ln Ai) while the entropy converges to a number of 
the same order at late times. Indeed, for Ai 1, the late-time entropy of entanglement has approximately the value 
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(34) 



which is dominated by the InAi term (the Ao term dies out at late times as long as Ap. is finite). This indicates that 
the constant Ai is due to the entanglement between the detector and the infinitely many degrees of freedom of the 
field. 

For Ai 3> 1, the effective temperature is approximately 



a 

of 



+ 7^(1 + 



7- 



$ 1 + 



7 + ifl 



1/2 



0(A- 



-l/2x 



(35) 



which, as we can see, is totally different from the Unruh temperature: It is determined not only by the proper 
acceleration a, but also by the properties of the detector (tt r ), the interaction (7), and the frequency cutoff of this 
theory (Ai). Moreover T g is very large (0(s/Ai)) for large Ai and non-vanishing even when a — > 0. 



B. Inertial detector: a — 



Substituting (|A1 1|) and (|A12|) into (l30|) . (I3T|) and (|33|) . one can write down in closed form the purity, entropy of 
entanglement and the effective temperature of an inertial detector (a = 0). These quantities are still nontrivial. In 
particular, in the ultra- weak coupling limit 7A1 <^ 1, one has 
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(36) 
(37) 



[Note that lim e „ >0 + ln(— 1 + e) = — Mr).] Both quantities go to zero as 7 
system returns to the nearly free theory. 



0, when the description of the combined 



C. Conventional result: the Markovian regime 



As we mentioned in Section III. A, the transition probability calculated from TDPT over infinite time is valid only in 
the ultra- high acceleration (temperature) limit or the ultra- weak coupling limit, both corresponding to the Markovian 
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regime for this model (with the quantum field as an Ohmic bath) . In these limits the entropy of entanglement and 
effective temperature do have the approximate values expected in the conventional picture: 

(a) Ultra-high acceleration (temperature) limit, a 3> 7Ai,f2. In this limit, the late-time entropy of entanglement 
reads 



In- 



+ 1 + 0(0, 



and the effective temperature 



T, 



s ~Tu + ( Ai + 7, 



In + 0(0 



(38) 



(39) 



is very close to the Unruh temperature TJy = ha/2Trkg. 

(b) Ultra-weak coupling limit, 7A1 <C a, n. In this limit, the late-time entropy of entanglement and effective 
temperature read 
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(40) 
(41) 



Again the effective temperature of the detector is very close to the Unruh temperature. 

It is well-known that the quantum mechanical time-dependent perturbation theory invokes the equivalent of a 
Markovian approximation (e.g., as in the derivation of the Pauli master equation [131]). Indeed, both the above 
limits correspond to the Markovian regime of the quantum Brownian motion, its dynamics being described by the 
master equation of Caldeira and Leggett [bij (in contradistinction to the non-Markovian master equation of Hu, Paz, 
and Zhang [||). There is admittedly a fluctuation-dissipation relation between the fluctuations in the field and the 
dynamics of the detector and thus backreaction is present, but with ultra- weak coupling it is restricted to the linear 
response regime (which has limited domain of applicability compared to the full noncquilibrium dynamics of the 
combined system). Here, the thermal bath is only slightly affected by the back reaction from the detector to the field, 
namely, the detector acts essentially as a test particle in the field. It is only under these special assumptions that the 
detector can come to equilibrium with an approximate thermal field at the Unruh temperature. Beyond this regime, 
the dynamics of the detector-field system has a totally different character. 



V. REMARKS 

We conclude with a few remarks: 

a. Working range of Unruh-DeWitt detectors. Conventional results from time-dependent perturbation theory 
over an infinitely long duration of interaction are trustworthy only in the ultra-high acceleration (or temperature) 
limit and the ultra-weak coupling limit, under the Markovian regime. In the non-Markovian regimes, the large 
constants Ai and Ao alter the scenario fundamentally. With the presence of the large constant Ai, which reflects 
the entanglement between the detector and the infinitely many degrees of freedom of the field, the evolution of the 
reduced density matrix of the detector collapses (sudden rapid change), as shown here, rather than diffuses (smooth 
gradual change), as conjured in the conventional picture. Furthermore, the detector at late times never sees an exact 
Boltzmann distribution over the energy eigenstates of the free detector, let alone the Unruh temperature. The strong 
interplay between the detector and the field makes the late-time quantum states for the detector-field system highly 
entangled and never factorizable (an analogy in quantum optics is the photon-atom bound state [15j|.) 

In equilibrium conditions, one can diagonalize the reduced density matrix and define an effective temperature (|35[) . 
which is usually large (~ 0(\/A~i~)) and does not vanish even for an inertial detector in Minkowski vacuum as long as 
the coupling is on. Only in the Markovian regime indicated above could this effective temperature get very close to 
the Unruh temperature. 

b. Correspondence with QBM. For a uniformly accelerated UD detector in (3+l)D with proper acceleration a 
the Unruh effect [H, 0, Q attests that it should behave the same way as an inertial UD detector in contact with a 
thermal bath at Unruh temperature Tjj, or more precisely, as an inertial harmonic oscillator in contact with an Ohmic 
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bath at Tu Q ■ This is clear from an examination of the integral in the derivation of the two-point functions of the 
detector, for example, Eq.(60) in Ref. [J]: 

Therefore for an inertial UD detector in contact with an Ohmic bath at temperature T, one can simply substitute 
2irkBT/h for a in the above results to get its purity, entropy and effective temperature. 

However, examining this from the vantage point of the exact solutions we obtained, we see the above statements 
are accurate only at the initial moment. After the coupling is switched on, the quantum state of the field will have 
been changed by the detector, so the field is no longer in the Minkowski vacuum and it does not make exact sense to 
say that the detector is immersed in a thermal state (or any state defined in the test-field description, i.e., where the 
field is assumed not to be modified by the presence of the detector) . 

A theorem by Bisognano and Wichmann (BW) [l6| states that the Minkowski vacuum, which is uniquely charac- 
terized by its invariance under all Poincare translations, is a Kubo-Martin-Schwinger (KMS) state with respect to 
all observables confined to a Rindler wedge. One may wonder why it does not apply here. The reason is that the 
BW theorem refers to the vacuum state of a quantum field alone, not the combined detector-field system. Even when 
the combined system is in a steady state, the quantum state of the interacting field is not invariant under spatial 
translations in Minkowski space, hence is not covered under the assumption of the BW theorem pertaining to Poincare 
invariance. 

Actually the Planck factor in (|42]) is a consequence of the BW theorem. Nevertheless, it is derived from only the 
free-field-solution part of the complete interacting field (see Eqs. (28), (56) and (58) in Ref. [4|). Here the factor is 
not distorted by the interaction simply because the field is linear and the coupling is bilinear. For nonlinear fields or 
couplings it would have a nonPlanckian spectrum and the departure from the conventional picture would be more 
pronounced. 

c. Backreaction and memory effects. It is common knowledge in nonequilibrium statistical mechanics [17l | that 
for two interacting subsystems the two ordinary differential equations governing each subsystem can be written as 
an integro-differential equation governing one such subsystem, thus rendering its dynamics non-Markovian, with the 
memory of the other subsystem's dynamics registered in the nonlocal kernels (which are responsible for the appearance 
of dissipation and noise should the other subsystem possesses a much greater number of degrees of freedom and are 
coarse-grained in some way). Thus inclusion of back-action self-consistently in general engenders non-Markovian 
dynamics. For our problem the two subsystems are the detector and the quantum field. Combining Eqs. (28), (30) 
and (10) in Ref.Q, we can write down the equation of motion for the detector evolution functions, 



(a»+n2) ff <+>(r;k) = ^. 

m 



^ +) (z(r);k) + A / dr'G Ict (z(r); z(r'))q^ (r' ; k) 



(43) 



which is an integro-differential equation. The backreaction to the detector is registered through the retarded Green's 
function G ro t of the field. Various approximations are usually invoked to solve this equation, amongst which the 
most common is the Markovian or memoryless approximation. This is what enters in the conventional derivation of 
the Unruh effect, but as we see above, is a very special and nongeneric subcase. 

d. Hiding Ai ? One may wonder why the large constant Ai cannot be absorbed by any parameter of this theory 
so one can renormalize T e g. To begin with, as we mentioned in Sec lIIlB) the UD detector theory is not a fundamental 
theory to meet the renormalizability requirement, and the presence of cut-offs as physical parameters is an expected 
feature which characterizes the range of validity of this semiclassical theory, just like the Compton wavelength of 
the electron acting as a cut-off in quantum optics. Second, the large constant Ai is not present in the renormalized 
stress-energy tensor of the field induced by the detector Q , thus T e g may not be a directly measurable quantity. The 
interference between the vacuum fluctuations and their back reaction cancels this cut-off dependence so that Ai is 
not observable outside of the detector. 

In some cases, though, Ai term can be subtracted from the physical quantities of the detector after these quantities 
are worked out. For example, one can subtract Ai from the detector energy defined in Eq.(82) of Ref. [H because 
during a physical process only the difference of energy matters. 

e. Hawking effect. It is tempting to see if analogous descriptions can be made and implications drawn for the 
Hawking effect. For a UD detector at rest in a static gravitational field, the response of the detector is similar. The 
simplest example is the UD detector fixed at radius r far from a Schwarzschild black hole, while the scalar field 
is in Hartle-Hawking (HH) vacuum (which is the counterpart of Minkowski vacuum for the uniformly accelerated 
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detector in Minkowski space :3:]). The response function per unit time for a massless scalar field in the Schwarzschild 
background is given by Candelas (Section V in Ref.[l9[). According to [l9[, when r — > oo, the response function 
for HH vacuum is exactly the same as the one in the uniform acceleration case (Eq.(58) in 0) with the Unruh 
temperature a/2irkB replaced by the Hawking temperature Th = l/SirksM for the Schwarzschild black hole with 
mass M . Therefore the results in this paper can be directly applied to the case of UD detector fixed at r — > oo outside 
a Schwarzschild black hole by neglecting its back reaction to spacetime. The effective temperature T e ff read off from 
such a detector in steady state is given by (|3"3")1 with a replaced by 1/4M rather than the Hawking temperature Th- 
For the case with the detector sitting very close to the event horizon, due to the effective potential barrier in the 
radial equation for the field, one has to consider the back-scattering of the retarded field induced by the detector in 
addition to the vacuum fluctuations described by the response function. We expect that the effective temperature 
read off from the detector in steady state would be even more complicated but definitely different from T H . We hope 
to report on this investigation in a later paper. 
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APPENDIX A: TWO-POINT FUNCTIONS OF INERTIAL DETECTORS IN MINKOWSKI VACUUM 



Recall that for a uniformly accelerated UD detector with proper acceleration a once we choose the factorized 
initial state (J2]), the two-point functions split into two part, (...) = (... ) a + ( ... )„. The coincidence limits of the 
two-point functions of the detector with respect to its ground state (with qo in Ref. [J] being zero) read 



™^ e- 2 ' m [fi 2 - 7 2 cos 2n V + jil sin 2fir?] , 
6{r,)e- 2 ^ \n 2 . - 7 2 cos 2^ - 7 ft sin 2^1 , 



(Al) 
(A2) 



which are independent of a, while the coincidence limits of the two-point functions of the detector with respect to 
Minkowski vacuum read 



( Q(v) 2 



= lim ~<{Q(r?),Q(r/)}K 

»7'->?j I 



2^7 
irmoQ' 



-%)Re { (a - In ^) e" 2 ^ sin 2 Qrj 



+ 



F 7+i0 (e-^) 
7 + in + a 



in 

7 



7 — 



(e- Q,) ) 



7 + iQ — a 



1 



< Q(v) 2 



^ + e -*m ^ + 1 - e'^Jj wrcoth^(ft - i 7 ) 

mno n^ Re {( 



(A3) 



2fvy 



Ai - In - Hr + A - In - e 



(n cos nr] — 7 sin nrj) 



+ ^(7 + zft) 2 e-^+ Q )" 



+ 4(7 +m 2 



in 



in 

1 

27?) 



7 + in — a 



7 + in + a \ 7 

'in 

1 



-2-rq 



1 + e 



-2iQ,r/ 



1 + e- 2iQr > i7rcoth-(f7 - i 7 ) 



(A4) 



Here again 77 = r — tq is the duration of interaction, 7 = Ag/87rmo, n = -y/fi 2 — 7 2 , F s (y) is defined by the hyper- 
geometric function as 



F s (y) = 2 F 1 (l + l,l,2+ -;y 
V a a 



(A5) 
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and 



4> s 



(A6) 



is the digamma function. The large constant Ao = — In £7 1 To — t' q \ — j e with the Euler's constant j e corresponds to 
the time scale of switching-on the interaction, so Ao could be finite in real processes, and for every finite value of Ao, 
the terms containing Ao in (|A3[) and (|A4[) vanish as 777 — » 00. The other large constant Ai = — lnO|r — t'\ — j e 
corresponds to the time-resolution or the cut-off frequency of this theory. Note that here we use slightly different 
definitions of Ai and Ao from those of Ref.[J| to make the arguments in logarithm functions dimensionless. Note also 
that the above results are valid when 0|r — t'\ <C 1 and H,\tq — Tq| <C 1. Beyond this regime, the form of the above 
"coincidence" limit of two-points functions should be modified. 

At late times (777 3> 1), the two-point functions of the detector with respect to the initial ground state ((. . .) a ) die 
away, and the two-point functions of the detector saturate to 



P 2 



Q 2 



Q 



27rmoil 
firriQ 



Re 



1 a 



2 



2iip. 



7+ifJ 



Re HO - 17) 



7 + id 



2iip, 



7+ifJ 



-Hm j (Ai - ln^- 



(A7) 
(A8) 



which arc identical to the results of the quantum Brownian motion of a harmonic oscillator in contact with an Ohmic 
bath at the Unruh temperature initially Q. To satisfy the uncertainty principle for all a (see FIGfT]), here we keep 
the constant Ai , which was subtracted in Ref. [1] because of the observation that Ai will not be seen outside of the 
detector. 

When a — » 0, the regular terms in the two-point functions diverge as In a and cancel the ones following Ao and Ai, 
so the two-point functions remain well-behaved. In this limit the two-point function (IA3|) continuously approaches 



Q(v) 2 



\a=0 



in 



6(r])Re {2Aoe- 2lri sin 2 

in 



7 



r [0, (7 + in) v ] - e - 2 ^ 



in 

7 



,-2jri 



in 

7 



1 



7 

a -2ifirj 



+ 1 _ e -2in„ )r[o,-(7 + *n)»/] 

n 



In 



7 



sin 2fi?7 



(A9) 



where T is the incomplete gamma function, while (| A4[) goes to 

/)/) 



QivY 



\a=0 



TrmoO 2 

in 



9(ri)Re ^ 



A ± n 2 



7 



r [0, (7 + in)n] + e 



-2-y-q 



A e 

in 



- 2lv (ncosn77-7sinn?7) 2 

l + e- 2in 'i ) r[0,-(7 + in)ry] 



in 

7 

-277) 



+ e 



-2717 



in 

1 



7 

-2ifirj 



In 



{i + inf 



n 



L 7 



-{n z + 7 2 ) + 27ncos2n?7 + (n 2 - 7 2 ) sin 2n?7 



At late times, one has 



Q 2 

p2 



hi 



\a=0 



■In 



7 



in 



27TTOon 7 + in 



2n 



(n 2 - 7 2 )ln 



7' 



in 



7 + in 



+ 7 



2Ai - In 1 



7_ 

n 2 



(A10) 

(All) 
(A12) 



Substituting this into (|30[) , (|3ip and (|33[) , one obtains the purity, entropy of entanglement and effective temperature 
for the incrtial detector in the Minkowski vacuum. 
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